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ABSTRACT 



Context. The Ara OB la association is a nearby complex in the fourth Galactic quadrant where a number of young/embedded star 
clusters are projected close to more evolved, intermediate age clusters. It is also rich in interstellar matter, and contains evidence of 
^{ , the interplay between massive stars and their surrounding medium, such as the rim Hll region NGC 6188. 

. Aims. We provide robust estimates of the fundamental parameters (age and distance) of the two most prominent stellar clusters, 

' NGC 6167 and NGC 6193, that may be used as a basis for studing the star formation history of the region. 

Methods. The study is based on a photometric optical survey (UBV1H„) of NGC 6167 and NGC 6193 and their nearby field, 
complemented with public data from 2MASS-VVV, UCAC3, and IRAC-Spitzer in this region. 
Results. We produce a uniform photometric catalogue and estimate more robustly the fundamental parameters of NGC 6167 and 
Q ' NGC 6193, in addition to the IRAS 16375-4854 source. As a consequence, all of them are located at approximately the same distance 

» I ' from the Sun in the Sagittarius-Carina Galactic arm. However, the ages we estimate differ widely: NGC 6167 is found to be an 

, intermediate-age cluster (20-30 Myr), NGC 6193 a very young one (1-5 Myr) with PMS, H„ emitters and class II objects, and the 

^ ■ IRAS 16375-4854 source is the youngest of the three containing several YSOs. 

Conclusions. These results support a picture in which Ara OB la is a region where star formation has proceeded for several tens of 
Myr until the present. The difference in the ages of the different stellar groups can be interpreted as a consequence of a triggered star 
formation process. In the specific case of NGC 6193, we find evidence of possible non-coeval star formation. 
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1 . Introduction getic activity of the stars in NGC 6167 and a giant HI bubble sur- 



. , , , , • , rounding was proposed by Arnal et al. (1987) and Waldhausen 

^ Ara OBla has been suggested as an example of triggered star ^^^^^^ ^^^^ p. ^ ^^j,^ 2008b). 

formation in the local Galaxy. A summary of the history of in- „, ■ ■• j- , • , xt^^ ^i 

„,. . J ,. J . The ages and distances of the two maior clusters NGC 6167 

vestigations of this region, our present understanding, and possi- j xT^^^ir.-, j ji r , , , j i i 

„° ^ ,. . ° \ J • xir 11 .1 /^rvrvo 1- j c and NGC 6193 dcrivcd thus far have been based On the photo- 

ble future studies is presented in WoUc et al. (2008a,b and refer- , j, i i r r , , 

^ electric and photographic photometry of a few stars, rather than 
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ences therein). This nearby association contains three open clus- , ^^t^ , , ^ cm w-l^v^v^^ . j ^^t^ 

,, „i 1 r , ■ ■ 1 1 on modern CCD data. Carraro & Munari (2004) presented CCD 

, ters (NGC 6193 in the center of the association, together with , „ xt^^ ^^r,^ i i - i ■ i j - i i- r 

" ^Tz-^^^i^T J xT/-^^ ^.-irv/ix T jj-.- 1 1 jj J 1 data for NGC 6204, but this cluster IS located in the outskirts of 

rrt NGC 6167 and NGC 6204). In addition, several embedded clus- , . . , i i j i • 

.^y. , • 1 , 1 \c ■ J • the association and may not be related to the main star/cluster 

ters are present together with both star-forming and quiescent j- • t^- • xt^^ ^irv-, r 

, , , r- 1 • rtr „ 1 ^r,r,o jT- m - i- f ormation proccss. Distancc cstimates to NGC 6 1 93 range from 

clouds (see Fig 1 in Wolk et al 2008a and Fig ^ m this pa- ^ ^ ^ ^^^^ .^^ ^ 3 ^^^^^^^ 

per). The idea that triggered star fromation has taken place in ^ q^,. p;',,^,,,^ , q^^^ T^,it.h»,™ 1 n.^.„^J^ qqi- x/,',.„», % 
Ara OBla is suggested by a variety of observations. The promi- 
nent cluster NGC 6193 appears to be connected to RCW 108 
(=NGC 6188), a rim HII region located westward that marks the 
edge of a molecular cloud being eroded by the hottest stars in the 



1973; Fitzgerald 1987; Kaltcheva & Georgiev 1992; Vazquez & 
Feinstein 1992). As for NGC 6167, the estimated distance ranges 
from 600 to 1200 pc and the age from 10 Myr to 40 Myr (Molfat 
& Vogt 1975; Bruck & Smyth 1967; Rizzo & Bajaja 1994). 



cluster, where several IRAS sources have been found (Comeron , " P'^'^f P^""' P^^^^^^ ^ CCD-based opti- 
et al. 2005). Star formation in the molecular cloud may have f ^ (UBVIc and //„) survey in an attempt to provide more ro- 
been triggered by NGC 6193. A connection between the ener- ^^'^.'^""f *^ and distance of NGC 6193 and 
NGC 6167. The photometry we present is sufhciently deep to 

* Member of Carrera del Investigador CONICET, Argentina allow us to estimate the nuclear and conti'action age, by detect- 

** On leave from Dipartimento di Astronomia, Universita di Padova, ing pre-main sequence (PMS) stars with //„ emission. Therefore 

Vicolo Osservatorio 2, 1-35122, Padova, Italy we expect to provide firmer estimates of the age and at the same 
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Fig. 1. Second generation Digitized Sky Survey (DSS-2; red filter) images, a) Image centered at ayaooo = 16 : 38 : 00; 6j2qoq ~ 
-49 : 15 : 20 and covering a 2° x 2° field. Solid rectangles centered in NGC 6167 and NGC 6193 indicates areas surveyed with 
UBVI observations. Dotted circle around NGC 6167 and several dotted circles around NGC 6193 indicates the areas with Ha(on/off) 
data, b) and c) Detailed images centered respectively at NGC 6193 and NGC 6167 showing main objects in each region and the 
adopted comparison fields. For the meaning of the symbols see in advance Fig|4]caption. 



time highlight possible differences between nuclear and contrac- 
tion ages. 

The layout of the paper is as follows. In Sect. |2] we describe 
the data sources, and their reduction and calibration. In Sect. |3] 
we present the analysis of the data together with the derivation 
of the main clusters parameters. In Sect. we derive the main 
properties of NGC 6167, NGC 6193, and the IRAS 16375-4854 
source. Finally, in Sect.|5]we discuss and summarize our results. 

2. Data 

2.1. Optical data 

Our main data set is based on CCD photometric observations of 
stars in the region of NGC 6167 and NGC 6193 (see Fig. □ car- 
ried out during several observational runs, and is complemented 
with available information from the literature. These data in- 
clude: 

- Las Campanas Observatory (LCO) UBVIc images of the re- 
gion of NGC 6167 acquired during two runs (25 June 2006 
and 24 June 2009) using the 1.0 m Swope telescope equipped 
with the Site 3 2048 x 3150 CCD camera. The field of view 
(FOV) of these images is about 14^8 x 22'8 with a scale 
a.'435/pix. The typical FWHM of the data was about a'9 and 
airmas values during the observation of the scientific frames 
were in the range 1.07-1.12. 

- Cerro Tololo Inter- American Observatory (CTIO) UBVIc 
images of the region covering about 40'0 x 40'.Q around 



NGC 6193 acquired during two runs (two frames on 25 and 
26 March 2006 and two frames on 21 May 2006). We used 
the Y4KCAM camera attached to the 1.0m telescope, which 
is operated by the SMARTS consortiurrQ. This camera is 
equipped with an STA 4064 x 4064 CCD with IS/J pixels. 
This set-up provides direct imaging over a FOV lO'.O x 20'0 
with a scale of 0'.'289/pi?(B Typical FWHM of the data was 
about 0'.'9, and airmasses during the observation of the sci- 
entific frames were in the range 1.02-1.28 . 

- Complejo Astronomico El Leoncito (CASLEO) Ha-on 
(656.6 nm) and Ha-off (666.6 nm) observations. They were 
performed on 14, 15, 17 and 18 May 2002 in six circular 
frames (4'5 radius) covering the area of NGC 6167 (1 frame) 
and NGC 6193 (5 frames). We used a nitrogen-cooled de- 
tector Tek-1024 CCD and focal reducer attached to the 215- 
cm telescope (scale O'.'813/pix). Exposure times in Ha-on and 
Ha-off ranged from 2 i to 150 i. 

- Other data sources: photometric data for a few bright stars 
that were found to be saturated in our observations, were 
taken from Webda databas^E 



Details of all these observations are given in Table [T] 



' http ://w w w. as tro .y ale. edu/smarts 

^ http://www.astronomy.ohio-state.edu/Y4KCarn/detector.html 
' http://www.univie.ac.at/webda 
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Table 1. Journal of observations of the scientific frames together 
with the derived photometric cahbration coefficients. 



Exposure times [sec] 


Observ. 


Frames U B 


V Ic N 


LCO 


long - 900 


600 600 1 




short 4x30 2x20 


10 10 1 


CTIO 


long 1500 1200 


900 700 4 




medium 200 100 


100 100 4 




short 30 30 


30 30 4 




vshort 10 7 


5 5 4 


Calibration and extinction coefficients 


Coefficient 


LCO 


CTIO 


Ml 


+4.528 ± 0.006 


+3.254 ± 0.007 


M2 


+0.125 ±0.011 


-0.034 ± 0.013 


M3 


+0.490 ± 0.010 


+0.450 


bi 


+2.938 ± 0.003 


+2.015+0.021 


b2 


+0.050 ± 0.004 


+0.120 + 0.025 


h 


+0.250 ± 0.010 


+0.250 


I'l 


+2.864 ± 0.003 


+1.802 + 0.013 




-0.068 ± 0.004 


-0.048 ± 0.016 


V3 


+0.160 + 0.010 


+0.160 


h 


+3.757 ± 0.007 


+2.676 ± 0.020 


h 


+0.039 ± 0.009 


-0.019 + 0.021 


h 


+0.080 + 0.010 


+0.080 







NGC 6167 


1 i 




K < 13 5 
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V<17 
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b) 
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2 4 6 
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Fig. 2. Radial density profiles for NGC 6167 (upper panel) and 
NGC 6193 (lower panel). Lowerplots (blue) correspond to CCD 
data (y < 16), whereas upper plots (red) correspond to 2MASS 
data {K < 13). Poisson error bars (Icr) are also shown. 



Notes. A' indicates the number of fields observed for each cluster region. 



Table 2. Approximate V magnitud limit values for 80% and 50% 
completeness factors for each studied band. 



Filter 


NGC 6167 


NGC 6193 




80% 


50% 


80% 


50% 


U 


13 


14 


16 


17 


B 


20 


21 


18 


19 


V 


21 


23 


21 


23 


I 


21 


23 


21 


23 


Ha 


13 


15 


13 


15 


J 


16 


18 


17 


19 


H 


16 


18 


17 


19 


K 


16 


18 


17 


19 


[3.6] 


12 


15 


13 


16 


[4.5] 


12 


15 


13 


16 


[5.8] 


12 


15 


13 


16 


[8.0] 


12 


15 


13 


16 



field. Finally, all data from different filters and exposures were 
combined and calibrated using DAOMASTER (Stetson 1992). 

2.1.2. Photometric calibration 

Standard stars were selected from particular areas of the catalog 
of Landolt (1992) (Mark A, PG 1323, and SA 110 for LCO 
data and SA 101, SA 104, and SA 107 for CTIO data) were 
used to determine the transformation equations relating our 
instrumental magnitudes to the standard UBVIq system. The 
selection of the fields was made to ensure that the stars had a 
wide range in colours. Aperture photometry was then carried 
out for all the standard stars (~ 70 per night) using the IRAF 
PHOTCAL package. To tie our observations to the standard 
system, we use transformation equations of the form 



2.1.1. Reduction 

All frames were pre-processed in a standard way using the 
IRAlQpack age CCDRED. For this purpose, zero exposures, and 
sky flats were taken every night. Photometry was performed us- 
ing IRAF DAOPHOT and PHOTCAL packages. Instrumental 
magnitudes were obtained using the point spread function (PSF) 
method (Stetson 1987). Since the FOV is large, a quadratic spa- 
tially variable PSF was adopted and its calibration for each im- 
age performed using several isolated, spatially well distributed, 
bright stars (about 25) across each field. The PSF photome- 
try was aperture-corrected for each filter and exposure time. 
Aperture corrections were computed by performing aperture 
photometry of a suitable number (about 20) of bright stars in the 

* IRAF is distributed by NOAO, which is operated by AURA under 
cooperative agreement with the NSF. 



u = U + uy+U2{U -B) + UiX (rms = 0.04) (1), 

b = B + bx+ b2(B -V) + b3X (rms = 0.03) (2), 

V = y -H v'l -H V2(B -V) + V3X (rms = 0.02) (3), 

i^Ic + ii + iiiV - Ic) + hX (rms = 0.02) (4), 

where UBVIc and ubvi are standard and instrumental magni- 
tudes respectively and X is the airmass of the observation. The 
used transformation coefficients and extinction coefficients are 
shown at the bottom of Table [T] 

The Ha-on/off observations were calibrated by fitting the em- 
pirical MS (Didelon 1982 and Feigelson 1983) to the adopted 
MS stars in each cluster field (see Sect. 13.41 and Fig. |7]l. In this 
way, the Ha(on) - Ha(off) index (hereafter A//„) has a nega- 
tive value for objects with a //„ line above the continuum. The 
relation between the AHa index and the //„ line width (Whu) is 
given by AH^ = -2.5 log(l - WHa[A]/60). 
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Fig. 3. Proper motion VPDs for adopted cluster areas of 
NGC 6167 and NGC 6193 and their corresponding comparison 
fields. Filled (blue) circles are adopted as likely members for 
each cluster 



2.2. Mid-infrared data 

The NGC 6193 and NGC 6167 regions were observed in the 
mid-IR by the Infrared Array Camera (IRAC) on-board the 
Spitzer Space Telescope (SST). IRAC was used to obtain im- 
ages in four channels (3.6, 4.5, 5,8 and 8.0 jdm) and a catalogue 
was produced using several of these photometric data points 
(GLIMPS^I). Observations covering NGC 6167 region were in- 
cluded in this catalogue. However, in the case of NGC 6193 and 
its surrounding area, only IRAC images were available in the 
SST database belonging to programs POOl 12, P00191, P03536, 
P20597, P30570, and P40321. 

We performed aperture photometry in all these images 
by using the IRAF PHOTCAL package following the proce- 
dure described by the Spitzer Science CenteiQ. Briefly, we run 
DAOFIND task to look for point sources. We then used the 
PHOT task using an aperture of 5 pixels and a background 
annulus from 12 to 20 pixels. Finally, to calibrate our results, 
we adopted the aperture correction values given in the on-line 
calibration tables. To avoid false detections, we selected only 
sources detected in both the four IRAC channels and the V band 
images. 

2.3. Compiementary data and astrometry 

Using the X-Y stellar positions obtained from our data, we 
correlated their positions with data from the following cata- 
logues: a) "Two-Micron All Sky Survey" (2MASS; Cutri et 
al. 2003, Skrutskie et al. 2006); b) "Vista Variables in the Via 
Lactea" catalogue performed at "Cambridge Astronomy Survey 



' http://www.astro.wisc.edu/glimpse/docs.html 
* http://ssc.spitzer.caltech.edu/irac/ 



Unit'CI (VVVh-CASU; Minniti et al. 2010), and c) "The Third 
U.S. Naval Observatory CCD Astrograph Catalog" (UCAC3; 
Zacharias et al. 2010). 

Our adopted procedure to perform the astrometric calibra- 
tion of our data was explained in Baume et al. (2009). The rms 
of the residuals in the positions were ~ 0'.'17, which is about 
the astrometric precision of the 2MASS catalogue (~ 0'.'12). As 
for NIR photometry, only stars detected in V filter were selected 
from NIR catalogues. We adopted 2MASS values only for the 
brightest stars (K < 13) and VVV-fCASU ones for the remain- 
der 

2.4. Data compieteness 

The completeness of the observed star counts is a relevant issue, 
thus it was computed by means of artificial-star experiments on 
our data (see Carraro et al. 2005). We created several artificial 
images by adding at random positions a total of 20 000 artifi- 
cial stars to our true images. These were distributed with a uni- 
form probability distribution of the same colour and luminosity 
as the real sample. To avoid overcrowding, in each experiment 
we added the equivalent to only 15% of the original number of 
stars. Since in all bands were considered only detections with a 

V filter counterpart, we performed those experiments on the long 
exposure image of V filter computing the completeness factor as 
the ratio of the number of artificial stars recovered to the num- 
ber of artificial stars added. The corresponding factors for the 
other bands were computed in a relative way as the ratio of the 
amount of detected stars in each band to the amount of detected 
stars in V filter corrected by completeness (weighted by diferent 
area coverage if necessary). 

The computed values of the completeness factor for differ- 
ent V magnitude bins are listed in Table |2] We note that for 
NGC 6167, only short exposures were obtained in the U filter, 
which explains the poor completeness of its catalogue in this 
band. Notwithstanding the completeness is sufficient to ensure a 
reliable estimate of the cluster parameters. 

2.5. Finai cataiogues 

The above procedure allowed us to build two astrometric, pho- 
tometric (13 bands) catalogues covering the regions of both 
NGC 6167 (33193 objects) and NGC 6193 (39338 objects). 
These two catalogues constitute the main observational database 
used in this study. A solid analysis of the behavior of the stellar 
energy distributions (SEDs) can be carried out with this tool, 
thus preventing possible degeneracies in the photometric dia- 
grams and allowing us to obtain more reliable results. Both cat- 
alogues are available only in electronic form at the Centre de 
Donnes astronomiques de Strasbourg (CDS). They includes X- 

Y positions, 2MASS identification (when available), equatorial 
coordinates (epoch 2000.0), and optical (UBVI and Ha-„„/off), 
NIR 2MASS/VVV (JHK), and MIR IRAC-SST (3.6, 4.5, 5,8 
and 8.0 fj.m) photometry. 

3. Data analysis 

3.1. dusters centers and sizes 

Radial stellar density profiles provide an objective tool to deter- 
mine the size of a star cluster They were computed, in a first at- 
tempt, starting from the estimated position of the cluster centers 

' http://casu.ast.cam.ac.uk/vistasp/vvv 
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Fig. 4. Optical TCDs of stars located in NGC 6167 (upper panel) and NGC 6193 (lower panel) regions. U - B ws. B - V diagrams: 
filled (blue) circles are adopted likely members for each cluster, open (red) circles are likely Ha emitter stars; big (orange) squares 
and filled (grey) triangles are respectively considered probable class I and class II objects (see Sect 13.51 and also the comments 
for each cluster); light dots are considered as field stars. The solid (blue) curve is the Schmidt-Kaler (1982) ZAMS, while dashed 
(green) lines are the same ZAMS, but shifted along the reddening line (red) by the adopted colour excesses indicated above them. 
They correspond to the adopted values for the cluster stars (see Table|4|. B - V vs. V-I diagrams: Symbols have the same meaning 
as in plots a) and c), whereas solid curves (blue) are intrinsic colours for luminosity class V and III from Cousins (1978ab). Dashed 
(red) arrows indicate the normal reddening path (Ry = 3.1). 



obtained by visual inspection of the second generation Digitized 
Sky Survey (DSS-2, red) plates, using our brightest catalogued 
stars {V < 16-17 depending on the cluster) and all the brightest 
2MASS catalogued objects (K < 13 - 13.5 depending on the 
cluster). We then repeated the procedure several times, slightly 
changing in each case the position of the centers and using differ- 
ent magnitude bins. We finally adopted these cases that produce 
the most uniform fall off in both bands. According to the vari- 
ation observed in different experiments, we estimated that the 
clusters centers have ~ 0'5 uncertainty in each coordinate. 

The cluster radial density profiles were computed in the 
usual way (see Baume et al. 2004), namely calculating the sur- 
face density in concentric, equal-area rings. Since we applied 
this method to both the optical (V filter) and the 2MASS in- 
frared data (K filter), this allowed us to detect the possible in- 
fluence of interstellar material. In both clusters, noticeable dif- 
ferences between the profiles for each band are detected. As a 
result, the sizes were estimated taking into account mainly the K 
filter, where absorption effects are minimized. 

The adopted centers and radii are indicated in Table H and 
Fig.[T] Centers are close to the ones given by Dias et al. (2002) or 
by the SI MB AD database, although the radius values are larger 
than those compiled by Dias et al. (2002), which are based only 
on visual inspection. The radial density profiles are shown in 
Fig. |2] Our data completely cover the clusters and sample also 
a significant portion of the surrounding field. This allows us to 



select different zones in each region and study them separately to 
understand the behavior of the stellar population in the directions 
of the clusters and their neighborhood (see Fig. [1). 



3.2. Proper motion analysis 

To obtain additional information to confirm and support the 
membership assignment, we used kinematic data provided by 
UCAC3 and Hipparcos catalogues (the later one was necessary 
for few bright stars with unacceptable large errors in UCAC3). 
Therefore, we used a sample of stars with UCACmag < 14 and 
proper motion errors smaller than IQmas/yr within the radius 
adopted for NGC 6167 and NGC 6193. To identify the members 
of the clusters we adopted the method developed by Cabrero- 
Cano & Alfaro (1985). This statistical algorithm fits the ob- 
served proper-motion distribution with the sum of two Gaussian 
distributions, one elliptical for the field and one circular for the 
cluster. Moreover, we took into account the observed errors in 
the proper motions for individual stars, which was suggested by 
Zhao & He (1990). The parameters of the two Gaussian distri- 
butions could then be determined from the maximum-likelihood 
technique. A summary of the application of these methods can 
be found in Fernandez Lajus et al. (2011). We note that the 
present study differs from that developed by Dias et al. (2006) 
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Fig. 5. Optical CMDs of stars located in NGC 6167 (upper panels) and NGC 6193 (lower panels) regions. Symbols are the same 
as in Fig.|4] Dashed (green) curves are the Schmidt-Kaler (1982) empirical ZAMS and the MS path from Cousins (1978ab). Solid 
curves are Marigo et al. (2008) isochones. All the reference curves are corrected by the adopted colour excess and apparent distance 
modulus for each cluster (see Table H). 



who used a total sample of 416 stars extracted from the UCAC2 
catalogue for the clusters under consideration. 

Table [3] shows the results obtained for NGC 6167 and 
NGC 6193. Finally, Fig. [3] presents the vector point diagram 
(VPD) for each cluster region. The adopted members (see 
Sect. [33) produce a clear clump separated from the field stars. 

3.3. Optical photometric diagrams 

To perform the final membership assignment of the stars to the 
clusters, we analyze both: a) the individual stellar positions in all 
the photometric diagrams (e.g. Baume et al. 2004, 2006, 2009); 
and b) the obtained membership probabilities values computed 
from the kinematic information (see Sect. I3.2l l. In the case of 
NGC 6193, Ha(on/off) data were also considered (see Sect. 13.4b . 
The optical two-colour diagrams (TCDs) and colour-magnitude 
Diagrams (CMDs) of NGC 6167 and NGC 6193 are shown in 
Figs. E] and |5l 

Individual positions of all stars down to V = 14 were simul- 
taneously examined in all the diagrams. At fainter magnitudes, 
UCAC3 does not provide reliable solutions. This makes contam- 



Table 3. Kinematic results from brightest stars (V < 14) in 
NGC 6167 and NGC 6193 regions 





Parameters 


NGC 6161 


NGC 6193 




f^a cos Sc 


-0.85 


0.54 


Cluster 


1^6, 


-2.61 


-2.91 




a 


1.40 


1.57 






-3.72 


-3.66 






-3.85 


-6.16 


Field 


£i 


8.34 


6.76 




£2 


8.03 


9.12 




P 


0.24 


0.28 



Note: 0"ffcosij.,Mi,) and (/j,vcos«, ) are the cluster and 
field proper motion centers in mas/yr; cr and (Xj, St) 
are the cluster and field proper motion dispersions in 
mas/yr; p is the obtained correlation coefficient. 



ination by field stars more severe, preventing an easy identifica- 
tion of faint cluster members using this method. 

The above selection of stars for both clusters allowed us to 
fit the blue edge of the adopted members with a properly red- 
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Fig. 6. Detailed CMD of stars located in NGC 6193 region. Symbols are the same as in Fig |4] Dashed (green) curves are the 
Schmidt-Kaler (1982) empirical ZAMS; solid (blue and red) curves are Marigo et al. (2008) and Siess et al. (2000) isochones for 
z = 0.02. All the reference curves are corrected by the adopted colour excess value and apparent distance modulus (see Table|4}. 
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Fig. 7. V vs. A Ha index of stars located in NGC 6167 region (left panel) and placed in NGC 6193 region and its surrounding (right 
panel). Symbols are the same as in Fig.|4] Dashed (green) curves are the empirical MS form using Didelon (1982) and Feigelson 
(1983) parameters. These curves are corrected by the adopted apparent distance moduli for each cluster (see Table|4}. Dotted vertical 
lines indicate the adopted limit to separate //„ emitting objects (see Sect. l3.4l i. 



dened Schmidt-Kaler (1982) zero-age main-sequence (ZAMS) 
in the U - B versus (vs.) Z? - V diagrams. We adopted the well- 
known relation Eu-b/Eb-v - 0.72 -i- 0.05 Eb-v and obtained the 
corresponding Eb-v value for each cluster The B - V vs. V - I 
diagram is routinely used as a diagnostic to detect deviations 
from the normal reddening law in a given direction. In this case 
the cluster diagrams suggest a normal law {R = Av/Eb-v) for 
NGC 6167 and a marginally different behavior for NGC 6193 
confirming previous results of WoUc et al. (2008b). 

We then simultaneously fitted a ZAMS or main sequence 
(MS) to all the CMDs to obtain cluster distance moduli, adopt- 
ing 7? = 3.1 for NGC 6167 and = 3.2 for NGC 6193. Finally, 
in Fig |6^ we compared the observed upper MS of NGC 6193 
with theoretical isochrones for solar metallicity, mass loss, and 
overshooting (Marigo et al. 2008). Some scatter (likely caused 
by the presence of binaries and rapid rotators) is present in our 
data and this prevents a unique isochrone solution. However, it 



is still possible to obtain an estimate of the nuclear age. All the 
resulting parameters are summarized in Table 2] 



3.4. Ha photometric diagrams 

One of the aims of this work was to detect stars with Ha emis- 
sion. For this purpose we used the AHa index described in 
Sect. 12.1.21 Similar indices were used in previous works (Adams 
et al. 1983; Sung et al. 1998; Baume et al. 2003). 

Figure |2] presents the V vs. AHa diagram for all stars with 
measured AHa index together with the corresponding ZAMS in 
this plane. Most stars follow approximately the ZAMS location 
indicating that they have no significant //„ emission allowing us 
to perform the calibration mentioned in sect. 12.1.21 Therefore 
positive AHa index values correspond to absorption features and 
negative values to emission ones. 
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Fig. 8. Infrared CCDs of stars located in NGC 6167 (upper panels) and NGC 6193 (lower panels) regions. Symbols are the same as 
in Fig.|4] Solid curves are the intrinsic Koornneef (1983) colours for MS stars. 



To detect stars with, on average, high Ha emission, we iden- 
tified likely //„ emission objects with the brighest stars {V < 
15.5) that had A//„ - 3eA//ff negative values. 

As for the distribution of H„ emitters, we note in Fig. [T] an 
evident concentration about T.Q westward of NGC 6193. This 
concentration has a diameter ~ 2' and is located close to the 
NGC 6188 rim (ay2000 ~ 16 : 40 : 30; 5y2000 ~ -48 : 44 : 20; 
see Fig. [T}?), just at the edge of the SFO 80 cloud (Sugitani et 
al. 1991; Sugitani & Ogura 1994). This group of H„ emitters 
closely corresponds to the MSX emission peaks at 8 fim found 
inside the cloud by Urquhart et al. (2009). All these elements 
seems to show different phases of stellar formation processes 
that are currently on-going. A particular spatial coincidence be- 
tween a class II object (see Sect. I3.5l l and an //„ emitter was 
found very near to GSC 08329-03510. This star is placed just 
at the border of the HII region NGC 6188 and very close to the 
center of the concentration of Ha emitters (see Sect. 13.4b . 

3.5. Infrared photometric diagrams 

Contracting stars emit an excess of radiation in the IR relative 
to their photospheric emission. This IR excess is considered to 
be produced by thermal emission from the circumstellar ma- 
terial. Therefore, this kind of objects are grouped into differ- 
ent "classes" according to their SEDs (Lada 1987; Andre et al. 
1993). Thus, NIR and MIR photometric observations allow one 
to detect and distinguish this kind of sources. 

We adopted the selection scheme presented in Gutermuth et 
al. (2008). This method is based on the location of the observed 
objects in the {K- [3.6])o vs. ([3.6] - [4.5])o plane using the stan- 
dard dwarf star colours from Bessell & Brett (1988) and Flaherty 
et al. (2007) colour excess ratios to obtain de-reddened colours. 



We then confirmed the selection by examining the locations in 
the J - H vs. H - K diagram (see Fig.|8]l. 

This analysis allowed us to recognize class I and class 
II/III candidates among our data. Their spatial distribution (see 
Fig. [TJ?) reveals that they are associated with IRAS sources 
and/or with denser molecular clouds reinforcing the hypothesis 
that they are indeed YSOs. 

4. Intrinsic cluster properties 

4.1. NGC 6167 

NGC 6167 appears to be an intermediate-age cluster with signif- 
icant differential reddening across its surface (see its TCDs and 
CMDs). This is consistent with the previous polarimetric results 
found by Waldhausen et al. (1999). 

The age of the cluster was computed by fitting the distribu- 
tion of the stars in U - B and B - V CMDs with isochrones. We 
paid particular attention to the brightest star in the CMDs, which 
is HD 149019 = CPD-49 9414 (AO la star from SIMBAD). It has 
a high membership probability value from the kinematic analysis 
(see Sect. l3.2l i. and is placed near the cluster center (see Fig.[T];). 
The spectrophotometric method also infers for this star a dis- 
tance moduli of (Vo - Mv)sp = 10.5, which is consistent with 
the adopted value for the cluster As a consequence, this star is 
adopted as a likely cluster member and the best-fit isochrone is 
that for 20-30 Myrs. In addition, only few stars in the field of 
this cluster exhibit a sizeable //„ emission (see Fig. |7^), bring- 
ing very support to this age value. 

The second bright star in the cluster region is HD 149065 
= CPD-49 9421 (B2 V star from SIMBAD). This star seems 
to be a foreground star according to its location in the photo- 
metric diagrams. In the U - B vs. B - V diagram, in particu- 
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Table 4. Parameters of the analyzed stellar groups 



Stellar group 


Center 


Radius 


Eb-v 


R 


V-My 


Av 


Vo~Mv 


Age[Myr] 




';i'2ooo 


52000 


['] 












Nuclear Contrac. 


NGC6167 


16:34:36.0 


-49:46:00.0 


5.5 


0.80 


3.1 


13.0 ±0.2 


~ 2.5 


10.5 ±0.2 


20-30 


NGC6193 


16:41:24.0 


-48:46:00.0 


9.0 


0.37 


3.2 


11.8 ±0.2 


~ 1.2 


10.7 ±0.2 


~ 1 -3-5 


IRAS 16375-4854 


16:41:18.9 


-49:00:34.0 


~ 2 


~ 2.2 


3.4 


~ 18 


~ 7.5 


~ 10.5 


< 1 



Note: Radii adopted values to select the sources for analysis in each region. 



lar, it shows a colour excess smaller than the cluster stars (0.24 
and 0.80, respectively). Its spectrophotometric distance moduli 
is (Vo - Mv)sp - 10.1, which is significantly smaller than the 
one adopted for NGC 6167 (10.5). These diff'erences might be 
caused by the Galactic absorption in the cluster direction (see 
Neckel & Klare, 1980), which is systematically greater by ~ 2 
magnitudes (decreasing the mean brightness from less than 1 
mag. to almost 3 mag.) at ~ 1 kpc from the Sun. On the other 
hand, HD 149065 is almost at the border of the adopted clus- 
ter region (see Fig. [1];), so we can assume that it is probably an 
interloper 

With an age of 20-30 Myr, it is unlikely that the sources con- 
sidered as class I/II objects in Fig.[8]are real YSOs. These objects 
appear to be young based only on the MIR diagram (Fig. |8};), 
although their location in the NIR diagram (Fig. |8^) should be 
different to adopt them as YSOs. 
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4.2. NGC 6193 

This cluster looks very young. Figure|6^ presents the reddening- 
corrected CMD obtained for the adopted likely cluster members 
in the upper MS. Differential reddening was removed using in- 
dividual colour excesses from Fig.|4j; and adopting R - 3.2. The 
isochrone fitting method (see Sect. 13. 3b then provides the nuclear 
age (~ 1 Myr). This result is compatible with the latest spectro- 
scopic analysis of its brightest cluster member (star HD 150136; 
03H-06V, Niemela & Gamen 2005). 

The lower MS of NGC 6193 shows quite a significant My 
scatter at any colour There are also several objects in the clus- 
ter field, which were those identified as Ha emitters and a few 
as class II objects (see Sect. 13.51 ). Figure [SJ? presents the loca- 
tion of these objects in a un-corrected CMD. These objects are 
all located above (~ 1-2 mag) the cluster ZAMS. Both effects, 
spread and location above the ZAMS, strengthen the idea that 
we are in the presence of a young cluster (cf. Fig. 3 in Preibisch 
& Zinnecker 1999). The observed spread might have intrinsic 
causes such as binarity/multiplicity, the random distribution of 
the orientation of the accretion discs around single and multiple 
systems, as well as external causes such as differential redden- 
ing, field star contamination, and photometric errors. 

The lower MS shape allowed us to derive an independent 
estimate of the cluster age (the contraction age) by comparing it 
with theoretical PMS isochrones. We achieved this by adopting 
the models computed by Siess et al. (2000). We note that it was 
impossible to find a unique isochrone that fitted the entire lower 
MS, but that we had to use a family of isochrone encompassing 
the stars loci in the CMD. Despite this, the mean location of the 
lower MS (see Baume et al. 2003) could be fitted with a group 
of isochrones that correspond to an age value marginally older 
than the nuclear age (~ 3-5 Myrs). 

We are aware that the accuracy of the contraction age might 
only be considered as an order of magnitude. This is because 
the large dispersion of objects in the lower MS may also reflect 



Fig.9.y vs. y-/ diagrams of stai-s located in IRAS 16375-4854 
region and its corresponding comparison field. Filled (blue) cir- 
cles are adopted members for each clump. Curves have the same 
meaning as in Fig|4] 



a non-coeval process of stellar formation. This is reinforced by 
a) the presence of objects of different nature in that area of the 
CMD such as stars near the ZAMS, //„ emitters, and class II ob- 
jects; and b) the existence of several young groups in the vicinity 
of NGC 6193 (see Wolk et al. 2008a,b and Fig. [It). 

Nonetheless, in the region close to NGC 6193 there are 
only a few class I or class II objects, despite our identification 
of several likely sources as //„ emitters. Low mass sources in 
NGC 6193 would then be somewhat older than those associated 
with IRAS sources and molecular clouds. In addition, the small 
number of class I/II sources would indicated a small, if any, age 
dispersion in the cluster itself (see Fig.|6];). 

4.3. Embedded clusters 

4.3.1. IRAS 16362-4845/RCW 108-IR 

This well known IRAS source contains the remarkable embed- 
ded and massive star formation region RCW 108-IR, which has 
been extensively studied in several wavelength ranges (from 
Radio to X ray; see Arnal et al. 2005, Comeron et al. 2005, 
Comeron & Schneider 2007, WoUc et al. 2008a and references 
therein). Our observations cover nearly the entire region of this 
object. However, sources in our catalogues were selected based 
on their detection in the V band and no further information to 
improve previous results could be achieved. 
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Fig, 10. NIR and MIR CCDs of stars located in IRAS 16375-4854 region. Symbols are the same as in FiglH Solid curves are the 
intrinsic Koornneef (1983) colours for MS stars. 



4.3.2. IRAS 16375-4854- IRAS 16379-4856 

These two IRAS sources are placed ~ 15' southward of 
NGC 6193 and were partially covered by our optical survey. 
IRAS 16379-4856 has an associated CO peak (Arnal et al. 2003) 
whereas IRAS 16375-4854 is located in a diffuse optical nebula 
that has several associated X ray sources (Wolk et al. 2008a). Our 
reduced IRAC SST data indicate that both IRAS sources do not 
seem to be independent since the objects detected in these bands 
form a ring structure of about 10' in diameter, which is possibly a 

bubble centered at ay2000 ~ 16 : 41 : 37; ^72000 49 : 02 : 32. 

This structure may have been generated by some energetic event 
at that position. 

By inspecting the stars detected on both regions, we note that 
IRAS 16375-4854 harbors objects brighter than IRAS 16379- 
4856. By using the corresponding optical photometric diagrams 
(CMDs and TCDs; see Figs. l9landlT0k). we estimated its main 
parameters (see Table |4]i. They were obtained by selecting the 
most probable members checking simultaneously their individ- 
ual position on all the photometric diagrams and also the finding 
chart. Apparently this region is affected by a non-ISM interstel- 
lar reddening {R ~ 3.4) and has a distance similar to NGC 6193 
and RCW 108-IR. 

As for the infrared photometric diagrams of IRAS 16375- 
4854, the position of the selected sources as YSOs appears to be 
unlikely. They might be class II/III objects, although we need to 
analyze X-ray observations to confirm this assumption. 

5. Discussion and conclusions 

Previous investigations by our group (Vazquez et al. 1996; 
Baume et al. 2003) presented optical observations (typically 
UBVRI) of young star clusters and star forming regions. To re- 
duce the uncertainties, we complement these UBVRI data with 
Ha photometry, kinematic information, and near/mid infrared 
data. This yields more reliable values of reddening, distance, 
and age for NGC 6167, NGC 6193, and the IRAS 16375-4854 
source. As a consequence, we can establish that all these groups 
are placed at the same distance from the Sun (~ 1300 pc) in the 
Sagittarius-Carina Galactic arm, but that their ages differ widely 
(see Table 4). We find that: 

- NGC 6167 is found to be an intermediate-age cluster (~ 20- 
30 Myr); 

- NGC 6193 is very young (-1-2 Myr) with PMS objects at 
different stages. Ha emitters, and class II objects; 



- Finally, IRAS 16375-4854 is the youngest of the three and 
contains several YSOs. 

In the specific case of NGC 6193, PMS stars appear to be 
slightly older than MS stars, i.e. the contraction age is older than 
the nuclear age, which might be an indication of non-coevality. 

The main result of this work is the homogeneity of distances, 
and the differences between the ages of the Ara OB la clusters, 
which may place on firmer ground future studies of the interac- 
tions between the different stellar and interstellar components of 
the region, such as those suggested by Arnal et al. (1987). The 
results presented here support the picture in which Ara OB la is 
a region where star formation has proceeded for several tens of 
millions of years up to the present. 
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